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Abstract. Line-formation calculations in the Rayleigh- Jeans tail of the 
spectral energy distribution are complicated by an amplification of non- 
LTE effects. For hot stars this can make quantitative modelling of 
spectral lines in the near-lR challenging. An introduction to the mod- 
, elling problems is given and several examples in the context of near-lR 

Oh' line formation for hydrogen and helium are discussed. 

o 

lyj 1 Introduction 

Near-infrared (near-IR) spectroscopy of early-type stars is a relatively new field. 
y—^ • It is motivated by the desire to study early phases of massive star formation 
^ (when the stars are deeply embedded in their parental clouds) and to investigate 
. young stellar populations throughout the Galactic disk and at the Galactic centre 
in particular. Despite hot stars showing a steep decline in flux levels towards 
' longer wavelengths, observations in the near-IR become highly competitive with 
optical region spectroscopy for regions of high extinction. In some cases near- 
IR observations are the only means to penetrate the dust. The most prominent 
' example is the Galactic centre, where the extinction in the if-band amounts to 
3 mag, compared with more than 30 mag in V. 

Near-IR spectroscopy became feasible only after substantial developments in 
detector technology over the past 15 years. Systematic observational studies of 
(normal) hot stars in the near-IR are nonetheless scarce. A spectral catalogue at 
^ . low resolution (i? = A/AA < 3000) was compiled by Hanson et al. (|1996p for OB- 
type stars. Classification spectra of early-type stars in the J, H and K bands 
can also be found in the atlases of Wallace & Hinkle (|1997p . Meyer et al. (|1998p 
and Wallace et al. pOQOp . More recently, intermediate- resolution spectroscopy 
(i?« 10 000) became the standard (FuUerton & Naiarro lTMSl Hanson et al. [2005)1 . 

Unguided by observations, the field attracted little interest on the theoretical 
side, except for an early prediction of photospheric Bra emission in B-typc stars 
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by Auer & Mihalas (|1969all . Later, Zaal et al. (|1999p could explain the observed 
emission cores of many near-IR hydrogen lines in early B-type stars via non-LTE 
effects, though several discrepancies between models and observation remained. 
The systematic behaviour of near-IR hydrogen and helium lines with spectral type 
and luminosity were investigated on theoretical grounds by Lenorzer et al. (|2004p 
for O stars. Finally, the most comprehensive study of near-IR spectra of OB 
stars and a comparison with analyses in the optical was performed by Repolust et 

al. mm . 

Considerable work has also been done on the massive star population near the 
Galactic centre. For the most part, this comprises extreme stars like Luminous 
Blue Variables and Wolf-Rayet stars (e.g. Najarro et al. 119971 120091 Martins et 
al. l2007p . However, such objects are beyond the scope of the present discussion as 
their line spectra are formed in the stellar wind. 



2 Challenges of non-LTE line-formation in the near-IR 



The challenges of near-IR line formation in hot stars can be summarised in brief. 
Let us recall the expression for the line source function. 



(fo,/6„)exp(Wfcr)-l 



(2.1) 



where 6/ and 6„ are the departure coefficients of the lower and upper level of the 
transition, respectively, v is the transition frequency, T is the temperature and h, 
c and k are the Planck constant, the speed of light and the Boltzmann constant, 
respectively. The line source function is identical to the Planck function, _B^, in 
LTE, i.e. deep in the stellar atmosphere {bi = = 1) or in cases where the lower 
and upper levels are tightly coupled (6/ = 6„). Spectral absorption lines experience 
non-LTE strengthening whenever Si^/ < \ (for 6/>6„), or they are weakened 
in the opposite case, potentially turning into emission lines when S'l exceeds the 
Planck function in the continuum. 

The source function is particularly sensitive to variations in the ratio of the 
departure coefficients 



|A5l 



bi/bu - exjp{-hiy/kT) 



{bi/bu - 1) + hv/kT 



(2.2) 



when hv/kT is small, i.e. in the Rayleigh- Jeans tail of the energy distribution. The 
denominator on the right hand side of Eqn. l2.2l can adopt values close to zero, am- 
plifying effects of a varying A(6//6„) considerably ('non-LTE amplification'). This 
can make near-IR lines in hot stars very susceptible to even small changes in the 
atomic data and details of the calculation. A quantitative reproduction of observed 
lines requires not only accurate knowledge of the plasma parameters (temperature 
and particle densities) and the radiation field in the stellar atmosphere, but also 
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Fig. 1. Observed hydrogen Unes (thick sohd lines) horn three series in Deneb (A2Ia): 
(a) Balmer, (b) Paschen, and (c) Pfund. Several strong metal lines are also present. 
The Paschen series is contaminated by many weak telluric features. Synthetic profiles 
from two model atmospheres are shown for comparison: a wind model (dashed line) 
and a hydrostatic model (dotted line). Note the discrepancies between observation and 
models in particular for the Pfund lines. From Aufdenberg et al. I|2002[l . reproduced by 
permission of the A AS. 

account for all relevant processes in model atoms and use of high-precision atomic 
data - even more so than required for analyses of hot star spectra in the visual 
region. 

3 Applications 

3. 1 Hydrogen 

The scope of problems with line-formation calculations in the near-IR can be exem- 
plified vifith one of the brightest stars of the northern hemisphere, the A-supergiant 
prototype a Cyg (Deneb) . Aufdenberg et al. (|2002p encountered problems with the 
simultaneous non-LTE modelling of the Balmer and near-IR hydrogen lines (see 
Fig. [I}, concluding, 'These failures indicate that a spherically symmetric, expand- 
ing, steady state, line-blanketed, radiative equilibrium structure is not consistent 
with the conditions under which [. . . ] the higher Pfund lines form.' In the following 
we want to investigate whether a less drastic explanation can be found. 
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Fig. 2. Comparison of effective collision strengths in hydrogen for several transitions n- 
n' , as indicated. The curves are: Butler (as published in Przybilla & Butler [20041 solid 
line), Johnson (fT972l dotted), Mihalas et al. (|1975l dashed), Percival & Richards (fT978l 
dash-dotted). From PrzybiUa & Butler (f2004ll . 
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Table 1. Model Atom Implementations for Hi 



Model Electron-impact excitation data 



A 


MHA (all n,n') 


B 


J72 (all n,n') 


C 


ABBS (n,n' <5), MHA (rest) 


D 


ABBS (n,n'<5), J72 (rest) 


E 


PB04 (n,n' < 7), PR (n,n' > 5), MHA (rest) 


F 


PB04 (n,n' < 7), PR (n,n' > 5), J72 (rest) 



ABBS: Anderson et al. H2000II : J72: Johnson H1972|l : MHA: Mihalas et al. H1975| l: PB04: Przybilla 
& Butler l(2004}; PR: Percival & Richards l|1978|l 

As explained above, small details do matter for line-formation calculations in 
the near-IR. The H i model atom used by Aufdenberg et al. looks reliable at first 
glance, given that 30 levels and over 400 transitions are accounted for explicitly 
in the non-LTE modelling. Moreover, the quantum-mechanical problem of the 
hydrogen atom can be solved analytically, such that the atomic data are of highest 
precision. But are they indeed? 

The previous statement applies to radiative data. Hydrogen is exceptional as 
one important source of systematic uncertainty compromising studies of most ions 
of the other elements is thus excluded. However, collisions by electrons intro- 
duce a third particle to the Coulomb problem, such that ab-initio calculations are 
required to determine reliable cross-sections for coUisional processes (excitation, 
ionization). Calculations of comprehensive sets of collisional data for all atomic 
levels of relevance in stellar atmospheres (lines up to n ~ 30 may be observed in 
supergiants, see e.g. Fig.[Tl while the Inglis- Teller limit is reached around n ~ 15 to 
20 for stars close to the main sequence) are still beyond present-day capabilities. 
However, it may be surprising to learn that reliable data for transitions involving 
levels of intermediate n were unavailable for a long time. 

In the absence of reliable data one has to resort to approximations. Approxi- 
mation formulae are based on simplifying theoretical considerations and indirect 
experimental evidence (see e.g. discussion by Mihalas [1967^ . As a boundary con- 
dition they should reproduce experimental constraints in the few cases where such 
are available. The approximations should provide data accurate to a factor better 
than two, and some studies claim the uncertainties to be as small as ~20%. 

A comparison of various approximation formulae with ab-initio data for effec- 
tive collision strengths for electron-impact excitation of several transitions having 
initial/final quantum numbers n, n' in neutral hydrogen is shown in Fig. [2] Good 
agreement is found between the frequently-used approximations according to John- 
son (|1972p and Mihalas et al. (|1975p with recent ab-initio data for the Lya transi- 
tion (1-2), as expected. The semi-empirical results of Percival & Richards (|1978p 
are displayed in this case for completeness only, as they are valid for n, n' > 5. 
The comparison for transitions with higher n, n' shows that uncertainty estimates 
for some approximations may be too optimistic. 

The question of which data should be preferred for non-LTE model atoms has 
to be answered by comparison with observation. Przybilla & Butler ([2004)) con- 
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Fig. 3. Spectrum synthesis for Br7 in Vega (AO V): models A, C, E and F (dotted, dash- 
dot-dot-dotted, full, dashed lines) and an LTE profile (dashed-dotted) are compared with 
observation (thick full line). Models B and D are omitted for clarity as they are almost 
identical to model F. From Przybilla & Butler (|2004p . 



structed several H i model atoms for this task (see Table [T] for a selection of data 
used for the evaluation of collisional excitation rates). Coverage of a wide range 
of plasma parameters for the tests is essential for the comparison, therefore, main 
sequence and supergiant stars of spectral types A, B and O were investigated. 
Stellar parameters were adopted from earlier detailed studies of their optical spec- 
tra. Note that the model atoms from Table [1] give almost indistinguishable profiles 
for the lines of the Balmer series. A few examples from the comparison arc now 
provided. 

The AOV star Vega is one of the most intensely studied stars, and in most 
respects well described by the assumption of LTE. The Balmer and Paschen lines 
are, in general, well matched by LTE computations, and these are practically 
identical to non-LTE results. However, comparison with the observed Br7 profile 
in Fig. [3] shows that LTE modelling fails to reproduce the line core, predicting a 
too shallow core. Non-LTE computations can improve on this, except for model 
A, which differs only slightly from LTE. Model A can therefore be disregarded for 
all other comparisons. However, on the basis of this one case alone a decision as 
to which of the model atoms B through F should be favoured cannot be drawn, 
as despite being noticeable the differences between these models are not highly 
significant. 

Non-LTE effects become stronger in supergiants, such that differences in the 
predictions from different model atoms can be expected to be amplified as well. 
Comparisons of predictions for BrlO and Brll (which are of photospheric origin. 
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Fig. 4. Higher Brackett lines in 77 Leo (AO lb) and /3 Ori (B8Ia): models E (full lines), 
F (dashed) and LTE calculations (dashed-dotted) are compared with observation (thick 
full). Profiles from models B and D resemble the model F predictions and are not shown 
here. Model C gives a slightly less good fit than model E. From Przybilla & Butler ([2004}. 



unaffected by the stellar wind) in the supergiants 77 Leo (AO lb) and (3 Ori (B8 la) 
are shown in Fig. H) Indeed, all models using the Johnson's (|1972p approximation 
formula (B, D and F) can be discarded as they produce too deep line cores. Equiv- 
alent widths are predicted too strong, by about a factor 2 in the more luminous 
supergiant /3 Ori (LTE equivalent widths are too low by about the same factor) . 
The remaining two models, C and E, produce similar predictions, with model E 
providing a slightly better fit to the observation, thus it becomes the recommended 
model for further use. The whole comparison process was more complicated than 
sketched here, involving more stars and more lines of the Balmer, Paschen, Brack- 
ett and Pfund series. Full details can be found in Przybilla & Butler (2004). 

The mechanisms driving departures of hydrogen from detailed balance in the 
atmospheres of early-type stars seemed to be well understood since the seminal 
work of Auer & Mihalas (1969b, 1969c), and numerous subsequent contributions 
- for line formation in the IR e.g. by Zaal et al. (|1999p . Indeed, various choices of 
the (mostly approximate) collisional data produce no significant differences in the 
stellar continuum or the Balmer line profiles, i.e. the features that are the starting 
point for quantitative analyses using model atmosphere techniques. On the other 
hand, the local processes that modify the radiatively induced departures from LTE 
become a dominant factor for line-formation computations in the IR. 

The example of /3 Ori is used to deepen the discussion. Departure coefficients 
bi=ni/n* (the and n* being the non-LTE and LTE populations of level i, 
respectively) for selected levels in calculations using models A-F are displayed in 



8 



Title : will be set by the publisher 



0.40 




Pj 
Er-/ 
I "7 



0.30 



3- 1C : 

4- 1G - 



0.20 



O 



0.1 ? 



0.00 



0.1 



-4 



3 



2 

log T, 







ross 



Fig. 5. Run of departure coefficients bi in /? Ori as a function of Rosseland optical depth, 
Tross, models A through F (dotted, dashed-dotted, dash-dot-dot-dotted, long dashed, full, 
dashed lines). The individual sets of graphs are labelled according to the levels' principal 
quantum number; all graphs for Hll coincide. Line-formation depths (core to wing) for 
a few features are indicated. From Przybilla & Butler (|2004[) . 

Fig. [5] The overall behaviour, i.e. the over- and under-population of the levels of 
the minor ionic species and the major species Hii, is governed by the radiative 
processes, while the differences in the collisional data lead to modulations. These 
are very small for the ground state and become only slightly more pronounced 
for the n = 2 level, as these are separated by comparatively large energy gaps 
from the rest of the term structure. Only colliding particles in the high-energy 
tail of the Maxwell distribution are able to overcome these energy differences at 
the temperatures in the stellar atmosphere. Thus, computations of the model 
atmosphere structure will not be significantly influenced, as the important bound- 
free opacities of hydrogen vary only in a negligible way. The IR line formation will 
be affected, as maximum effects from variations of the collisional data are found 
for the levels with intermediate n at line-formation depth. 

In fact, the differences in the collisional cross-sections from approximation for- 
mulae and ab-initio computations are largest for transitions among the n — 3-7 
levels with An — 1 and 2, and they can amount to more than an order of magni- 
tude. The higher Rydberg states show less sensitivity as they approach the limit- 
ing case of LTE, which is independent of the details of individual (de-)populating 
mechanisms. Collisional cross-sections from ab-initio computations up to n ~ 7 
are therefore sufficient to eliminate a significant source of systematic error. Using 
the available data it turns out that the MHA- (models A, C, E) and J72-type 
approximations (models B, D, F) give rise to basically two sets of distinct be- 
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Fig. 6. Ratio of line source function Sl to Planck function at line centre as a function 
of Tross in /3 Ori. Line designations as in Fig. O From Przybilla & Butler (|2004p . 



haviour, with the former tending to dampen non-LTE departures more efficiently 
than the latter, due to larger collisional cross-sections. Such differences in the 
level populations are the cause for the line-profile variations as in Figs. [3] and ID 
In fact, the bi typically vary by only several percent, but non-LTE amplification 
changes equivalent widths by a factor ~2 in the example shown in Fig. |4]. The 
corresponding reaction of the line-source function to variations in the collisional 
data is displayed in Fig. [51 

Finally, we may ask how the comparison of spectrum synthesis based on the 
recommended model atom E looks like for Deneb, the starting point of our dis- 
cussion. Indeed, an excellent match of model and observation can be obtained, as 
shown by Schiller & Przybilla pOOSp . see Fig.[7l Note that the lower members of 
the Balmer, Paschen and Brackett series are affected by the stellar wind of this su- 
pergiant. Consequently, computations based on the hydrodynamic non-LTE code 
Fastwind are shown for the comparison in this case, while hydrostatic calcula- 
tions with Detail/Surface based on an Atlas9 atmosphere are sufficient for 
the lines of photospheric origin. 

3.2 Helium 

As in the case of hydrogen, non-LTE effects on He i lines are stronger in the near- 
IR than in the optical. An example for the (relatively weak) A2.11 /im singlet and 
triplet feature in the if-band spectrum of r Sco (B0.2 V) is shown in Fig.[8l These 
lines experience a very pronounced non-LTE strengthening. The stellar parameters 
for the calculations were determined by analysis of the optical spectrum. 

Details of the model atom and of the modelling assumptions also play an im- 
portant role for non-LTE line-formation computations for He i in the near-IR. Con- 
siderable differences exist between photoionization cross-sections used in classical 
studies of Hei line formation in early B-type stars and modern ab-initio data. 
A comparison of the data for the metastable 2s '^S state (the lower level of the 
He 1 10 830 A transition) is shown in Fig. [9j displaying results of Gingerich (|1964p 
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Fig. 7. Modelling (black) of the hydrogen lines in the visual and the near-IR spectrum 
of Deneb (gray curve). Overall, a good to excellent agreement is achieved, compare 
with Fig. [1] The synthetic spectra are calculated with a hybrid non-LTE technique 
using Detail/Surface (photospheric lines) or Fastwind (FW, wind-affected lines), as 
indicated. Except for the Pfund series, all panels show the same range in relative flux. 
Some lines, such as P/3 and If/3, are noticeably affected by the stellar wind. According 
to Schiller & Przybilla (p008|) . 
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Fig. 8. Modelling of the He I A2.11/^m singlet and triplet feature (if-band) in rSco 
(B0.2 V). Like in the case of hydrogen, near-IR transitions of He I experience, in general, 
stronger non-LTE effects than the lines in the visual. From Nieva & Przybilla (|2007p . 




Fig. 9. Comparison of the photoionization cross-sections for the He I 2s ''S state (Fernley 
et al. 119871 full line), as used in our reference He I model atom, and of Gingerich l|1964l 
dotted), as used in previous model atoms by Auer & Mihalas (I1973P and Husfeld et 
al. (|1989|l . The cross-sections differ by more than 50% at threshold. From Przy- 
billa (|2005)l . 

and of Fernley et aL (|1987p . as obtained within the Opacity Project. The modern 
data give a cross-section at threshold larger by more than 50%. Note that the 
presence of the resonance structure at short wavelengths is irrelevant for the prob- 
lem, as the stellar flux is negligible there, resulting in an insignificant contribution 
to the photoionization rate. 
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Fig. 10. Comparison of observed equivalent widths for He 1 10 830 A in B-type stars 
(Lennon & Dufton 119891 filled circles; Leone et al. 119951 open circles) with model pre- 
dictions: Auer & Mihalas (|1973l dotted line), Dufton & McKeith (|1980l dashed line), and 
our own calculations (full lines) for microturbulence ^ = (lower) and 8kms~^ (upper 
curve). The abscissa is the reddening-free [ci] index, a temperature indicator, spanning 
a range from early A to late O. From Przybilla (|2005p . 

The too low cross-sections were already corrected in calculations by Dufton & 
McKeith ()1980p . leading to predictions of strong emission for He i 10 830 A in early 
B-type stars in contrast to absorption predicted by Auer & Mihalas fll973p using 
the Gingerich (|1964p data. Surprisingly, observations by Lennon & Dufton (|1989p 
favoured the predictions made with the incorrect data, at least qualitatively, as the 
line was found to be in absorption in most cases (Fig. [TOl) . A satisfactory repro- 
duction of the observed trend was also not obtained later by Leone et al. (|1995p . 
Shortcomings in the model atmospheres were invoked in both cases. 

Again, a small detail of the calculations has a large effect on the predicted line 
profile. Good agreement between the observed trend and model calculations for the 
near- IR line Hei 10 830 A was finally obtained by Przybilla (1 2005^ . using modern 
atomic data and accounting for line-blocking effects in the non-LTE computations, 
see Fig. [101 It turns out that the traditionally analysed He i lines in the optical 
region are not very sensitive to whether line blocking is accounted for or not. 

The effects of using photoionization cross-sections from Gingerich (|19641 as 
implemented in the He I model atom of Husfeld et al. I1989P and of Fernley et 
al. (|1987[ as implemented in the reference model by Przybilla 12005^ on departure 
coefficients of the n = l and 2 levels are shown in Fig. [TT] Only the departure 
coefficients for the 2s "^S state show differences of as little as ~5% at line-formation 
depths, i.e. only the Hei 10 830 A transition is affected - 2p^P° and 2p^P° are the 
lower levels of the He i lines in the optical region. Note that as a consequence line- 
formation depths also change. Almost the entire effect can be accounted for by 
the 2s ^S photoionization cross-section. Practically all other differences between 
the two sets of data are negligible for the problem. 
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Fig. 11. Comparison of departure coefficients bi as a function of Tross for the n=l 
and 2 levels of Hel from computations using the Husfeld et al. H1989I dotted lines) 
and our reference model atom (full lines), and computations for this new model atom 
with the photoionization cross-section for the 2s state replaced by that used in the 
older one (dashed lines). Line-formation depths for Hel A10 830A are indicated. The 
computations are for a stellar atmospheric model with Tea = 30 000 K, log (? = 4.0 (cgs), 
^ = Okms~^, solar metallicity, and solar helium abundance. From Przybilla l|2005p . 
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Fig. 12. Ratio of S'l to Bi, at the centre of the 10 830.34 A fine-structure component 
of the transition for the same atmospheric parameters as in Fig. 1111 In addition to 
line designations from Fig. 1111 the dashed-dotted line shows results from a computation 
using the new model atom, but neglecting line blocking. In the inset the resulting line 
profiles are compared after convolution with a Gaussian of FWHM of 20kms~^. Use 
of the improved photoionization cross-sections turns the absorption into emission, while 
neglect of line blocking gives excess emission. From Przybilla (|2005[l . 
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Fig. 13. Comparison of the normalised J-band spectrum of the extreme helium star 
V652Her (histogram) in the region around He I A10 830 A with three model predictions: 
LTE (dashed line) and two non-LTE computations using the model atom of Husfeld et 
al. (|1989l dotted line) and the improved model atom of Przybilla (|2005l full line). In 
contrast to the other transitions, He I A10 830 A experiences strong non-LTE weakening. 
An unbroadened synthetic spectrum is displayed at the bottom in order to illustrate what 
can be expected from high-resolution observation. From Przybilla et al l|2005[) . 



The response of the line source function to photoionization cross sections is 
shown in Fig. [TH Non-LTE calculations based on the (incorrect) older photoion- 
ization data yield absorption lines that are shallower than in LTE. Computations 
based on modern cross-sections give weak emission, while neglect of line blocking 
results in excess emission. Finally, the predicted trends discussed in the literature 
are thus understood. 

He 1 10 830 A is not the only near-IR line of neutral helium that is sensitive to 
non-LTE amplification. A comparison of the observed spectrum of the extreme 
helium star V652 Her in the J-band with synthetic spectra for two different model 
atoms (and an LTE calculation) shows that all lines are affected (Przybilla et 
al. 120051) . Our reference model atom succeeds in reproducing all the lines simul- 
taneously. While He 1 10 830 A experiences non-LTE weakening, the other features 
are subject to non-LTE strengthening. 

4 Conclusions & Outlook 

We have seen that amplification of non-LTE effects is ubiquitous for near-IR lines 
in hot stars. This makes non-LTE line-formation calculations challenging, as every 
aspect of the modelling has to be considered properly. On the other hand, the high 
sensitivity offers the possibility to put tight constraints on the atomic data selected 
for the construction of non-LTE model atoms. Deficits in the data will show up 
in the comparison of the model predictions with observation. 
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The near-IR studies of hot stars so far comprise relatively strong lines of hydro- 
gen and helium, based on medium-resolution spectra. The new generation of high- 
resolution near-IR spectrographs like Crires on the VLT promise to revolutionise 
the field. In particular, a multitude of (weak) metal lines will become accessible for 
the first time (see e.g. Przybilla et al. 120091 Nieva et al. I2009p . The development 
of proper non-LTE models will turn quantitative near-IR spectroscopy of hot stars 
into a crucial tool for Galactic studies, and with the upcoming extremely large 
telescopes (diffraction-limited observations using adaptive optics will be feasible 
only at near-IR wavelengths) also for extragalactic stellar astronomy. 
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